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Abstract 

Background: ParaHox and Hox genes are thought to have evolved from a common ancestral ProtoHox cluster or 
from tandem duplication prior to the divergence of cnidarians and bilaterians. Similar to Hox clusters, chordate 
ParaHox genes including Gsx, Xlox, and Cdx, are clustered and their expression exhibits temporal and spatial 
colinearity. In non-chordate animals, however, studies on the genomic organization of ParaHox genes are limited to 
only a few animal taxa. Hemichordates, such as the Enteropneust acorn worms, have been used to gain insights 
into the origins of chordate characters. In this study, we investigated the genomic organization and expression of 
ParaHox genes in the indirect developing hemichordate acorn worm Ptychodera flava. 

Results: We found that P. flava contains an intact ParaHox cluster with a similar arrangement to that of chordates. 
The temporal expression order of the P. flava ParaHox genes is the same as that of the chordate ParaHox genes. 
During embryogenesis, the spatial expression pattern of PfCdx in the posterior endoderm represents a conserved 
feature similar to the expression of its orthologs in other animals. On the other hand, PfXIox and PfGsx show a novel 
expression pattern in the blastopore. Nevertheless, during metamorphosis, PfXIox and PfCdx are expressed in the 
endoderm in a spatially staggered pattern similar to the situation in chordates. 

Conclusions: Our study shows that P. flava ParaHox genes, despite forming an intact cluster, exhibit temporal 
colinearity but lose spatial colinearity during embryogenesis. During metamorphosis, partial spatial colinearity is 
retained in the transforming larva. These results strongly suggest that intact ParaHox gene clustering was retained 
in the deuterostome ancestor and is correlated with temporal colinearity. 
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Background 

The ANTP superclass homeobox genes including Hox 
and ParaHox encode transcription factors that play cru- 
cial roles in many aspects of development in bilaterian 
animals [1]. ParaHox and Hox genes are thought to 
have evolved from a putative common ancestral gene 
complex, the ProtoHox cluster, prior to the divergence 
of cnidarians and bilaterians [2-4]. A recent study fur- 
ther proposed that distinct Hox and ParaHox loci were 
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present in the last common ancestor of all animals [5]. 
On the other hand, analysis on the cnidarian Hox/ 
ParaHox genes suggested that the ParaHox cluster 
formed as a result of tandem duplication rather than 
cluster duplication [6]. One of the best-known properties 
of Hox genes is the spatial and temporal colinearity 
between their expression patterns and their positions 
within clusters on a chromosome [7]. Hox genes at the 
3' end are involved in patterning the anterior of the em- 
bryo, genes in the middle of the cluster pattern the cen- 
tral regions of the embryo, and the genes at the 5' end 
pattern the posterior of the embryo in a phenomenon 
called spatial colinearity [8]. Temporal colinearity refers 
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to a situation in which anterior Hox genes are expressed 
earlier and posterior genes are expressed later [9]. These 
observations have led to the hypothesis that the physical 
organization of Hox genes on the chromosome is im- 
portant for proper morphological differentiation along 
the anteroposterior axis [10]. 

ParaHox cluster was first identified in the amphioxus 
Branchiostoma floridae, in which three member genes, 
Gsx, Xlox, and Cdx, are linked in a genomic region with 
Gsx adjacent to Xlox in the same orientation, followed 
by Cdx on the opposite strand [11,12]. This cluster or- 
ganization is conserved in Xenopus, mouse, and human 
[12,13]. Studies on the genomic organization of ParaHox 
genes in non-chordate animals are limited to only a few 
animal taxa and the clustering seems labile. In pro- 
tostomes, for example, C. elegans has only one Cdx 
ortholog, pal-1, while Drosophila has Gsx (ind) and 
Cdx orthologs that are not linked on the chromosome 
[14-16]. In the annelid Platynereis dumerilii, Gsx and 
Xlox genes are linked together, whereas Cdx is located 
away from these two genes [17]. In deuterostomes, the 
ParaHox genes are not linked in teleost fishes [18], the 
ascidian Ciona intestinalis [19], and the sea urchin 
Strongylocentrotus purpuratus [20], although all three 
genes are found. Recently, an intact ParaHox cluster was 
found in the starfish Patiria miniata [21], suggesting 
that ParaHox clustering was maintained in the echino- 
derm ancestor and the sea urchins have modified the 
original arrangement. 

The maintenance of ParaHox clusters in chordate ge- 
nomes also led to discussions about their spatial and 
temporal colinearity in gene expression. The expression 
of amphioxus ParaHox genes exhibits both temporal and 
spatial colinearity [11,12,22]. AmphiCdx transcript is the 
first among the three ParaHox genes to be detected by 
in situ hybridization, being expressed initially around the 
blastopore during mid-gastrulation and later in the larval 
hindgut and posterior neural tube. After gastrulation, 
AmphiXlox expression is detected in the posterior 
archenteron and neuroectoderm. In the larval stage, 
AmphiXlox transcripts mark the midgut-hindgut bound- 
ary. AmphiGsx is the last gene to be expressed, and its 
expression is restricted to a few cells in the anterior 
neural tube. Therefore, similar to the spatial colinearity 
of Hox genes along the anteroposterior axis, amphioxus 
ParaHox genes also exhibit spatial colinearity with Gsx 
expressed in the anterior neural tube, Xlox and Cdx 
expressed in the middle and posterior endoderm and 
neural tube, respectively. The temporal colinearity of 
amphioxus ParaHox genes is also evident, but inverted 
with respect to the pattern in the Hox cluster: transcript 
of the posterior Cdx is detected first and the anterior 
Gsx is expressed last [22]. Studies on the expression of 
ParaHox genes in other deuterostome and protostome 



animals including sea urchin [20], ascidian [23-26], mouse 
[27-29], polychaete worms [17,30,31], and gastropod [32] 
have also shown similar expression domains: Gsx genes 
are mostly expressed solely in the central nervous system 
(CNS) with a rostral anterior limit; Xlox genes are ex- 
pressed both in the CNS and the central regions of devel- 
oping guts, such as the pancreas of vertebrates; Cdx genes 
are expressed in more posterior regions of the CNS and 
gut. Temporal colinearity, on the other hand, is reversed 
in S. purpuratus and lost in C. intestinalis [33]. 

Hemichordates are the sister group of the echinoderms, 
which together are referred to as the Ambulacraria that is 
closely related to chordates. Hemichordates retain bilat- 
eral symmetric body plan throughout their life and share 
several morphological similarities with both echino- 
derms and chordates. Therefore, hemichordates have been 
served as a model to gain insights into the origins of 
chordate or deuterostome characters [34,35]. Hemichor- 
dates consist of two major groups, the solitary enterop- 
neust acorn worms and the colonial pterobranchs. The 
phylogenetic relationships among pterobranchs and two 
of the enteropneust groups, direct-developing Harriman- 
nidae and indirect- developing Ptychoderidae, have been 
uncertain. Fossils of both pterobranchs and enteropneusts 
were present in Middle Cambrian strata [36] . Phylogenetic 
analyses on 18S rDNA have placed pterobranchs as a 
sister group to the Harrimaniidae [37,38], making ente- 
ropneusts paraphyletic. A recent study that compared 
microRNA repertoires among several hemichordate and 
echinoderm species has unambiguously supported the 
monophyly of enteropneust worms [39]. The indirect- 
developing Ptychoderidae is of particular interest be- 
cause their tornaria larvae share striking developmental 
and morphological similarities with echinoderm larvae 
[40-42]. Ptychodera flava from the Ptychoderidae has re- 
cently been shown to possess a 12-gene Hox cluster that 
is similar to the organization of chordate Hox clusters but 
with different posterior genes [43]. In this study, we in- 
vestigated the genomic organization and expression of 
ParaHox genes in P. flava. We found that P. flava contains 
an intact ParaHox cluster with a similar arrangement to 
that of chordates. We also showed that expression of these 
genes exhibit temporal colinearity but lose spatial co- 
linearity during embryogenesis. Nevertheless, partial 
spatial colinearity is retained in the transforming larva. 
These results strongly suggest that intact ParaHox gene 
clustering is correlated with temporal colinearity and was 
retained in the deuterostome ancestor. 

Results and discussion 

A previous study has identified four P. flava ParaHox 
genes, which include PfGsx, two Xlox (PfLoxl and 
PfLox2), and PfCdx [44]. Our PCR analysis revealed that 
PfGsx, PfLox2, and PfCdx are located on a single BAC 
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clone, PfBSllFlO. However, the presence of the PfLoxl 
gene in the genome is questionable because a BLAST 
search of the P. flava genome (using shotgun sequen- 
cing; unpublished data) did not identify PfLoxl, More- 
over, several attempts to amplify PfLoxl from genomic 
DNA failed, and the expression of PfLoxl could not be 
detected in any developmental stages analyzed by RT- 
PCR and in situ hybridization. Furthermore, another 
indirect-developing hemichordate, Balanoglossus simo- 
densis, contains a single Xlox gene (BsXlox) [45]. PfLoxl, 
PfLox2, and BsXlox proteins share identical homeo- 
domains, and the overall identity between PfLoxl and ei- 
ther PfLox2 or BsXlox is comparable (83% and 80%, 
respectively) (Additional file 1: Table SI). Thus, it is pos- 
sible that, similar to other invertebrate deuterostomes 
[11,20,21,45], P. flava contains only one single Xlox gene 
(PfLox2). This conclusion may indicate that PfLoxl is an 
ortholog gene from a related species collected in the 
same area as P. flava and that the clone was obtained 
from a contaminated library made from pooled animals 
[44]. We hence renamed PfLox2 as PfXlox to represent 
the single Xlox gene in P. flava. 

Genomic organization of P. flava ParaHox genes 

The sequence data of the BAC clone PfBSllFlO was as- 
sembled into a single continuous sequence of 169 kb in 
length (Figure 1A). PfGsx and PfXlox are organized in 
the same transcriptional orientation with a 47.4 kb in- 
tergenic region, while PfCdx is located 6.2 kb away from 
PfXlox in a head-to-head orientation. The total genomic 
sequence of the P. flava ParaHox cluster occupies a re- 
gion of 106.8 kb. No other protein coding genes were 
identified within the cluster by routine sequence sear- 
ches. This cluster organization is remarkably similar to 
the intact chordate ParaHox clusters found in am- 
phioxus, Xenopus, mouse, and human [13]. No intact 
ParaHox cluster has been found outside deuterostomes. 
To explore the genomic organization of ParaHox genes 
in protostomes, we searched several recently published 
spiralian genomes. In the marine polychaete Capitella 
teleta genome, Xlox and Cdx are on the same scaffold 
with three predicted genes between them and Gsx is lo- 
cated on a different scaffold [30,46]. Two Gsx, one Xlox, 
and one Cdx genes are found in the leech Helobdella ro- 
busta genome, although each is located on one of four 
different scaffolds [46]. In the limpet Lottia gigantea 
genome, the Gsx and Xlox genes are linked together but 
Cdx is separated on a different scaffold [46]. The Pacific 
oyster Crassostrea gigas genome contains three ParaHox 
genes, with Gsx and Xlox located on the same scaffold 
and Cdx on another [47]. In the pearl oyster Pinctada 
fucata genome, Gsx, Xlox, and Cdx are located on three 
different scaffolds [48]. Therefore, currently there is no 
evidence for the existence of intact ParaHox clusters in 



protostome genomes (Figure IB). Searching the draft 
genome of Saccoglossus kowalevskii, a Harrimaniidae 
acorn worm, revealed that while the SkXlox and SkCdx 
genes are next to each other in the same head-to-head 
orientation, SkGsx is on a different scaffold. It remains 
unknown whether the three S. kowalevskii ParaHox 
genes form a cluster in the genome. Nevertheless, the 
close linkage of PfGsx, PfXlox, and PfCdx in P. flava sug- 
gests that ParaHox genes were clustered in the deutero- 
stome ancestor. 

It has been suggested that the ParaHox cluster of the 
chordate ancestor is flanked by CHIC and PRHOXNB 
genes [12]. The 5' flanking region of PfGsx in PfBSllFlO 
is too short to confirm if CHIC is located there in 
P. flava. Using a BLAST search against the whole genome, 
we identified a genomic contig containing P. flava 
PRHOXNB gene but found that this contig was not 
located on PfBSllFlO, suggesting that the P. flava 
ParaHox cluster is not flanked by PRHOXNB. On the 
other hand, FLT1, the ortholog of which has been shown 
to be a neighbor of human and mouse ParaHox clusters 
[12], was found next to PfCdx (Figure 1A). Therefore, as 
in the single mammalian ParaHox cluster, genomic reor- 
ganizations have occurred around the P. flava ParaHox 
cluster but not within it. 

Temporal colinearity of P. flava ParaHox genes 

To investigate the temporal order of P. flava ParaHox 
gene expression, we performed quantitative PCR (QPCR) 
at different embryonic stages (Figure 2). None of the 
ParaHox genes were detected in the unfertilized egg as 
maternal transcripts. The PfCdx transcript was first de- 
tected at the blastula stage (16 hours post fertilization; 
hpf) and peaked at the late gastrula stage (43 hpf). PfXlox 
transcript was barely detected in the blastula and early 
gastrula (22 hpf) stage and the expression level increased 
at the late gastrula stage. The PfGsx transcript was de- 
tected at a low level a few hours after PfXlox was expres- 
sed, and the transcript could not be detected at the 
tornaria larva stage (65 hpf). These results suggested that 
the first ParaHox gene to be expressed in P. flava is PfCdx, 
followed by PfXlox, and finally by PfGsx. As observed in 
amphioxus, the activation of P. flava ParaHox genes re- 
tains temporal colinearity. 

Spatial expression patterns of P. flava ParaHox genes 

The spatial expression domains of the three ParaHox 
genes were determined by in situ hybridization (Figure 3). 
Similar to its ortholog AmphiCdx in amphioxus, PfCdx 
transcript was the first among the three ParaHox genes to 
be detected by in situ hybridization around the blastopore 
at the early gastrula stage (Figure 3C, white arrow). In the 
mid-gastrula stage, after the mesoderm forms at the tip of 
the archenteron, PfCdx was strongly expressed around the 
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Figure 1 Genomic organizations of ParaHox genes. (A) Phylogenetic tree depicting genomic organization of the P. flava ParaHox cluster (top) 
compared to the clusters of amphioxus {Branchiostoma floridae) and mouse {Mus musculus) on chromosome (Chr.) 5. Blue, green, red, purple, and 
yellow boxes are exons of Gsx (Gsh), Xlox (Ipfl), Cdx, PRH0XNB, and Fltl, respectively. The light blue line indicates the fully sequenced PfBS1 1 F10 
BAC clone. Bent arrows indicate transcriptional orientations. The bottom line is a zoomed out view showing the neighboring Fltl gene located 
near the mouse ParaHox cluster. (B) Evolution of the ParaHox cluster in bilateria. The phylogenetic tree represents the genomic organization of 
ParaHox genes in several bilateria, including protostomes (Pr) and deuterostomes (De). Protostomes can be further divided into lophotrochozoa 
(Lo) and ecdysozoa (Ed). The following lophotrochozoa animals were included in the analysis: two polychaete species {Platynereis dumerilii and 
Capitella teleta) and one leech species {Helobdella robusta) belonging to the Phylum Annelida (An); limpet {Lottia gigantea) and two oyster species 
{Crassostrea gigas and Pinctada fucata) in the Phylum Mollusca (Mo). Two ecdysozoa species fruit fly Drosophila melanogaster and nematode 
Caenorhabditis elegans do not contain the full complement of ParaHox genes. In deuterostomes, Phylum Echinodermata (Ec), including sea urchin 
Strongylocentrotus purpuratus and starfish Patiria miniata, and Phylum Hemichordata (He), including Ptychodera flava and Saccoglossus kowalevskii, 
constitute Ambulacraria (Am) that is closely related to Phylum Chordata (Ch). ParaHox gene organizations of the three chordate species, 
amphioxus Branchiostoma floridae, ascidian Ciona intestinalis, and mouse Mus musculus, are presented. IDs of the scaffolds (Sc) on which ParaHox 
genes are located are indicated beneath the illustrated scaffolds found in the genome databases. Blue, green, and red triangles indicate the 
orientations of Gsx, Xlox, and Cdx, respectively, in the genome. Double slashes between two genes indicate that although the genes are located 
on the same chromosome or scaffold they are separated by intervening genes. Species names shown in red contain intact ParaHox clusters. 



blastopore and the posterior archenteron (Figure 3D-E, 
white arrows). Later, at the tornaria larva stage, PfCdx ex- 
pression persisted in the ectoderm around the blastopore 
and the hindgut until at least 21 dpf (days post fer- 
tilization) (Figure 3F-H). The anterior boundary of the 
PfCdx expression domain demarcated the midgut-hindgut 
boundary (Figure 3F-H, black arrows). PfXIox expression 
was detected around the blastopore at the mid- and late 
gastrula stage (Figure 3L-M, black arrowheads). Expres- 
sion of PfXIox became weaker in the tornaria larva but 



remained in the same ectodermal domain surrounding the 
blastopore and the posterior endoderm (Figure 3N-P, 
black arrowheads). PfGsx expression was restricted to a 
few cells around the blastopore in the late gastrula stage 
and disappeared in the tornaria larva (Figure 3T-X, white 
arrowheads). To analyze spatial expression domains in de- 
tail, we performed double fluorescent in situ hybridization 
of PfXIox with either PfCdx or PfGsx (Figure 4). PfCdx 
transcript can be clearly observed around the blastopore 
at the mid-gastrula (MG) stage (Figure 4A-C). In the late 
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Figure 2 Temporal expression profiles of ParaHox genes. 

Transcript levels of ParaHox genes were measured by QPCR at 
different embryonic stages including unfertilized egg (0 hpf), early 
blastula (12 hpf), late blastula (16 hpf), early gastrula (24 hpf), late 
gastrula (43 hpf), and tornaria larva (65 hpf). The Y-axis is the 
transcript level normalized to the 18S rRNA level at the same stage. 



gastrula embryo (LG) and tornaria larva (TL), the ex- 
pression domains of PfCdx and PfXlox in the posterior 
ectoderm domain mostly overlapped (Figure 4D-I). PfCdx 
expression in the posterior endoderm extended to the 
midgut-hindgut boundary whereas PfXlox only expressed 
in the most posterior endoderm. The expression pattern 
of PfCdx in the developing gut conforms to the generally 
conserved feature as observed for its ortholog genes in 
other animals. On the other hand, PfXlox expression 
around the blastopore and the most posterior endoderm 
is different from the conserved expression pattern in the 
middle endoderm. The restricted expression of PfGsx also 
resided within the ectodermal domain around the blasto- 
pore with no detectable anterior staining (Figure 4J-L). 
These data indicated that expression of the P. flava 
ParaHox genes lose their spatial colinearity during the em- 
bryonic stages. The anterior boundary of the posterior 
PfCdx gene expression domain is relatively more anterior 
than that of the anterior PfGsx and middle PfXlox genes. 

The adult body plan of indirect developing acorn 
worms, such as P. flava, develops from a feeding tor- 
naria larva that undergoes metamorphosis [34]. We fur- 
ther investigated the expression pattern of P. flava 
ParaHox genes in the transforming larva collected from 
plankton tows. PfCdx expression was detected in the 
posterior endoderm (Figure 5A-C, arrows) in a pattern 



UF B EG MG LG TL TL TL 



0 hpf 16 hpf 22 hpf 30 hpf 43 hpf 3 dpf 9dpf 21 dpf 




Figure 3 Expression patterns of P. flava ParaHox genes. In situ hybridization analyses of PfCdx (A-H), PfXlox (l-P), and PfGsx (Q-X) were performed 
on unfertilized egg (UF), blastula (B), early gastrula (EG), mid-gastrula (MG), late gastrula (LG), and tornaria larva (TL) at different time points (hpf, hours 
post fertilization; dpf, days post fertilization). White and black arrows indicate the expression of PfCdx and the midgut-hindgut boundary, respectively. 
Black and white arrowheads denote the expression of PfXlox and PfGsx, respectively. The observed views are indicated at the bottom right corner of 
the panels in the first row. Unless otherwise indicated, the observed views of subsequent panels are the same as the first panel in each column. 
DV, dorsal view; LV, lateral view; W, ventral view. The scale bar is 50 urn. Panels T and U' are the blastopore view (BV) of panels T and U, respectively. 
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Figure 4 Double fluorescent in situ hybridization of PfXIox with PfCdx or PfGsx. DNP-labeled PfXIox probe and DIG-labeled PfCdx probe 
were used for in situ hybridization on mid-gastrula (A-C), late gastrula (D-F), and tornaria larva (G-l). The arrows and arrowheads indicate the 
anterior boundaries of the PfCdx and PfXIox expression domain, respectively. Double stainings of PfGsx and PfXIox transcripts were performed on 
the late gastrula embryos (J-L). 



similar to its expression at the tornaria larva stage 
(Figure 3). Surprisingly, during transformation of the 
competent larva, PfXIox was expressed in the middle 
part of the endoderm (Figure 5D-F, arrowheads). This 
expression pattern is similar to that of Xlox genes in 
other animals. The PfGsx transcript, however, could not 
be detected in the transforming larva. These data indi- 
cated that the expression patterns of PfXIox and PfCdx 
in the endoderm of the transforming larva conform to 
the spatial colinearity that is conserved with other ani- 
mals. Thus, we conclude that the spatial colinearity in 
the P. flava ParaHox cluster is altered during the em- 
bryonic stages, but remained in the transforming larva. 




Figure 5 Expression of PfCdx and PfXIox in transforming larvae. 

Whole mount (A, D) and sagittal sections (B, E) of the transforming 
larvae hybridized with antisense probes against PfCdx (A-C) or PfXIox 
(D-F). Panels C and F are higher magnifications of panel B and E, 
respectively. Arrows and arrowheads indicate the in situ hybridization 
signals of PfCdx and PfXIox, respectively. The scale bar is 200 urn. 



Evolution of the P. flava ParaHox cluster 

We have demonstrated that in the hemichordate 
P. flava, expression of the three ParaHox genes exhibits 
temporal colinearity during embryogenesis. However, 
despite having an intact ParaHox cluster, the spatially 
staggered expression pattern of these genes is lost in 
embryogenesis but partially maintained in the transfor- 
ming larva. These results reinforce the idea that regu- 
latory mechanisms controlling temporal colinearity are 
the major constraints that maintain an intact ParaHox 
cluster [7]. Spatial colinearity, on the other hand, seems 
to be independent of clustering. The Cdx expression pat- 
tern in the posterior endoderm is extremely conserved 
in bilateria regardless of whether Cdx is located within 
the cluster or not. Thus, the cis-regulatory module 
(CRM) that controls posterior endoderm expression of 
Cdx must be closely associated with Cdx genes during 
animal evolution. The embryonic expression of PfXIox 
and PfGsx, on the other hand, represent novel patterns. 
Both the conserved neuroectoderm and middle endoderm 
expression of Xlox genes are not detected in P. flava 
embryos; instead, embryonic expression of PfXIox is in the 
ectodermal cells around the blastopore. We hypothesize 
that this novel expression pattern is due to a new CRM 
generated during evolution. Compared to the genomic 
structures of the chordate Xlox genes, the relatively 
long intron between the two exons of the PfXIox gene 
provides the possibility for generating the novel CRM 
(Figure 1A). The conserved CRM controlling PfXIox 
expression in the midgut is still retained as this 
expression domain is maintained in the transforming 
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larva. Further studies on the cis-regulatory analysis of 
PfXlox will be required to test this hypothesis. 

A widely held view is that Hox genes are primarily 
used in larval cells destined to become parts of the adult 
body rather than in the development of larva-specific 
structures [33]. A similar idea, however, has not been 
discussed with respect to ParaHox genes. Our data show 
that the temporal colinearity of ParaHox genes is ob- 
served during embryogenesis whereas spatial features 
are more conserved in transforming larvae than during 
embryogenesis. These data seem to support the idea that 
ParaHox genes are used similarly to Hox genes in pat- 
terning anteroposterior axis for adult structures. Never- 
theless, the expression of sea urchin ParaHox genes 
exhibits spatial colinearity in embryonic gut despite ha- 
ving a broken cluster [20]. Further investigation on the 
roles of ParaHox genes in indirect-developers may help 
to solve this issue. 

Conclusions 

It is generally proposed that ParaHox and Hox genes 
have evolved from a common ancestral ProtoHox clus- 
ter. However, unlike the widespread existence of Hox 
clusters in the animal kingdom, intact ParaHox clusters 
have only been found in chordates and recently in a star- 
fish. Here we find that the hemichordate P. flava con- 
tains an intact ParaHox cluster without any intervening 
genes and the organization of the cluster matches that 
of chordates and the starfish. Our finding suggests that 
both the ancestral Ambulacraria and ancestral deutero- 
stome possessed intact ParaHox clusters (Figure IB). 
The current status of the protostome genome assemblies 
could not provide definitive evidence for the presence of 
a scaffold containing all three ParaHox genes. Improve- 
ment of the current genome assemblies or sequencing 
other protostome genomes will be required to reveal 
whether protostomes contain intact ParaHox clusters. 

The intact P. flava ParaHox cluster represents a spe- 
cial case in the discussion of temporal and spatial co- 
linearity. The expression of these genes exhibit temporal 
colinearity but spatial colinearity is modified. Our data 
suggest that gene clustering is correlated with temporal, 
not spatial, colinearity. These findings may encourage 
more intensive studies on the regulatory mechanisms for 
maintaining gene clusters. These results, together with 
the recent finding that Hox cluster of P. flava has similar 
organization to the chordate Hox clusters [43], reinforce 
the idea of using P. flava to gain insights into the origins 
of deuterostome and chordate body plans. 

Methods 

Construction of BAC libraries 

We constructed two sets of Ptychodera flava BAC libra- 
ries (PfBS and PffiH) estimated to represent the 4 (PfBS) 



or 5 (PfBH) x genome coverage for each library. BAC li- 
braries were constructed according to the procedures as 
described previously [49]. Briefly, genomic DNA was 
prepared from sperm of a single P. flava individual col- 
lected at the sand bar, Kaneohe bay, Oahu, Hawaii. Ge- 
nomic DNA was partially digested with Sac I (for PfBS) 
or Hind III (for PfBH) and subjected to pulsed-field gel 
electrophoresis. The DNA fragments corresponding to 
110 (for PfBS) or 135 (for PfBH) kb were isolated and li- 
gated into the pKS145 (for PfBS) or pKS200 (for PfBH) 
vector. Transformation was carried out by electropo- 
ration using E. coli DH10B as a host strain. Ampicillin- 
resistant transformants were collected and then stored 
in 384-format plates at -80°C 

Screening of BAC libraries 

To facilitate screening by PCR, BAC clones were orga- 
nized in three-dimensional pools using Biomek FX. 
DNA was prepared from each pool using Kurabo PI-200 
and PI-1100. Specific primers for each ParaHox gene 
were designed based on the homeobox sequences. Ini- 
tially, the pools of the PfBS library were screened by 
PCR using iyQfo-specific primers and two positive 
clones were identified. Subsequently, DNA prepared from 
the isolated BAC clones was used as a template, and PCR 
was performed to determine the presence of P. flava 
ParaHox genes in the clones. 

BAC DNA sequencing and assembly 

One P. flava BAC clone, PfBSllFlO, containing three 
ParaHox genes, PfGsx, PfXlox, and PfCdx, was sequen- 
ced using standard shotgun procedures at the National 
Institute of Genetics, Japan. Individual BAC DNA was 
isolated using Kurabo PI-200. End sequences of the 
shotgun library were analyzed in an Applied Biosystems 
3730x1 capillary sequencer, and reads were assembled 
using the Phrap-Consed suite programs. 

Identification of ParaHox genes in genome databases 

ParaHox gene searches were conducted using available 
genome databases. We used P. flava ParaHox sequences 
as queries to search the genomes of Capitella teleta (JGI 
vl.0), Helobdella robusta (JGI vl.0), Lottia gigantea (JGI 
vl.0), Crassostrea gigas (GigaDB), Pinctada fucata 
(Genome Ver 1.00), and Saccoglossus kowalevskii (NCBI 
Skowl.l). 

Animals, embryos, and larval collection 

Adult acorn worms {Ptychodera flava) were collected 
from the Penghu Islands, an archipelago off the western 
coast of Taiwan. Animals were induced to spawn using 
the temperature shift method as described [50]. Eggs 
were washed extensively with filtered seawater before 
fertilization [51]. Embryos were cultured at 23°C and fed 
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with Rhodomonas lens after they hatched from the 
fertilization envelope. Transforming larvae were col- 
lected by plankton tows at Sand Island offshore, 
Honolulu, Hawaii. Our experimental research was ap- 
proved by Academia Sinica Biosafety Review & Biomate- 
rials and Lab Biosafety Information System (certificate 
number BSF0410-00002036). 

Quantitative PCR (QPCR) 

Total RNA was extracted from several embryonic stages 
using Trizol reagent (Invitrogen). Ten micrograms of 
Trizol-extracted RNA was further purified using the 
RNeasy Micro Kit (Qiagen) to remove genomic DNA. 
One microgram of the purified RNA was reverse tran- 
scribed using the iScript cDNA synthesis kit (Bio-Rad). 
The resulting cDNA was used as a template for QPCR. 
Levels of 18S rRNA were used to normalize samples. 
QPCR analysis was performed on a Roche LightCycler 
480 with the LightCycler 480 SYBR Green I Master 
(Roche). The QPCR primers used in this study were 
designed based on the published sequences [44,52], and 
the primer sequences were as follows: PfGsx-QF: 5'- 
CGCTGTCAACAGTGCCTTAG -3'; PfGsx-QR: 5 - CA 
CCTTCGCAGAGAGAGGAA -3'; PfXlox-QF: 5 -GGC 
GG AC AACAAG AAC ACTT-3 PfXlox-QR: 5 -CCCGC 
AGAAAGACTGACTTC-3 PfCdx-QF: 5-ACTTTC 
GTCACGGCAGACTT-3 PfCdx-QR: 5-GTCCGAA 
TGGTGAGCTTGTT-3'; Pfl8S-QF: 5 -CCTGCGGCTT 
AATTTGACTC-3 Pfl8S-QR: 5 -AACTAAGAACGGC 
CATGCAC-3'. 

In situ hybridization 

cDNA clones of PfGsx and PfCdx were obtained from li- 
brary screening [44]. PCR cloning was used to construct 
a PJXlox cDNA clone from embryonic cDNA that was 
amplified with PfXlox-F (5 -CCAACATGGAGAGTTC 
TAATCC-3 ) and PfXlox-R (5 -GCGGTCTGTCTTT 
GTCAGAT-3') primers. Antisense riboprobes were syn- 
thesized from these cDNA clones using digoxigenin 
(DIG) RNA Labeling Mix (Roche) with T7 or SP6 RNA 
polymerases (Promega). For double fluorescent in situ 
hybridization, dinitrophenol (DNP) labeled probes were 
made using the LabellT DNP Labeling Kit (Mirus). Em- 
bryos were fixed and dehydrated as described [42]. After 
embryos were rehydrated in phosphate buffer saline 
containing 0.1% Tween-20 (PBST), they were digested 
with 10 ug/ml Proteinase K for 5 min, washed twice with 
0.2% glycine, and fixed with 4% paraformaldehyde in 
PBST for 1 hr. Embryos were then washed sequentially 
with 0.1 M triethanolamine, 0.25% acetic anhydride, and 
0.5% acetic anhydride. After washing extensively with 
PBST, hybridization was performed and the embryos 
were imaged as previously described [53]. For in situ 
hybridization of transforming larvae, similar procedures 



were applied except that the larvae were digested with 
10 (ig/ml Proteinase K for 20 to 30 min and the antibody 
was preabsorbed with adult powder [51]. The stained 
competent larvae were embedded in Tissue-Tek Optimal 
Cutting Temperature (O.C.T.) compound (Sakura, Japan) 
and cryosectioned at 10 um. 

Additional file 



Additional file 1: Table SI. Comparison of PfLoxl to other 
ambulacraria Xlox proteins. Accession number: Ptychodera flava Lox1 
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Balanoglossus simodensis Xlox (BsXIox), AB506760. Saccoglossus kowalevskii 
Xlox (SkXIox), XM_002741 106. Strongylocentrotus purpuratus Lox (SpLox), 
NM_2 14650. 



Abbreviations 

ANTP: Antennapedia; BAC: Bacterial artificial chromosome; BLAST: Basic local 
alignment search tool; CNS: Central nervous system; CRM: Cis-regulatory 
module; DIG: Digoxigenin; DNP: Dinitrophenol; dpf: Days post fertilization; 
hpf: Hours post fertilization; PBST: Phosphate buffer saline with Tween-20; 
PCR: Polymerase chain reaction; QPCR: Quantitative polymerase chain 
reaction; RT-PCR: Reverse transcription-polymerase chain reaction. 

Competing interest 

The authors declare that they have no competing interests. 
Authors' contributions 

YHS, MIA, JKY, HS, and NS conceived and supervised the project. Tl carried 
out sequence analysis. YCC performed gene expression analyses. RA 
performed double florescent in situ hybridizations on the Hawaiian acorn 
worm embryos. HCT and CHT participated in molecular cloning and larval 
culture, respectively. KT and TH collected fecund animals and transforming 
larvae from Hawaii. RK assembled the P. flava draft genome. AF established 
the P. flava BAC libraries. Tl and YHS wrote the manuscript. All authors read 
and approved the final manuscript. 

Acknowledgements 

The authors wish to thank Paola Oliveri from University College London and 
Kevin Peterson from Dartmouth College for providing the cDNA clones and 
Eric Rottinger for providing fixed P. flava embryos from the Hawaiian 
population for comparison with the population in Taiwan. This project was 
supported by KAKENHI on Grant-in-Aid for Scientific Research (C) from JSPS 
(23570266) and the Grant-in-Aid for Scientific Research on Innovative Area 
"Genome Science", from Ministry of Education, Culture, Sports, Science 
and Technology in Japan (to Tl), the National Science Council, Taiwan 
(101-2627-B-001-001 to YHS and 101-2923-B-001-004-MY2 to YHS and JKY), 
and a Career Development Award from Academia Sinica, Taiwan (to JKY). 

Author details 

Marine Genomics Unit, Okinawa Institute of Science and Technology 
Graduate University, Onna, Okinawa 904-0495, Japan, institute of Cellular 
and Organismic Biology, Academia Sinica, Nankang, Taipei 1 1529, Taiwan. 
3 Stazione Zoologica Anton Dohrn, Napoli 80121, Italy. 4 Marine Biological 
Laboratory, Graduate School of Science, Hiroshima University, Hiroshima 
722-0073, Japan. 5 Pacific Biomedical Research Center, University of Hawaii, 
Manoa, HI 96822, USA. 6 National Institute of Genetics, Yata 1111, Mishima, 
Shizuoka 41 1-8540, Japan, department of Biological Sciences, Graduate 
Schools of Science and Engineering, Tokyo Metropolitan University, 
Hachiohji, Tokyo 192-0397, Japan, institute of Biogeoscience, Japan Agency 
for Marine-Earth Science and Technology (JAMSTEC), 2-15 Natsushima-cho, 
Yokosuka 237-0061, Japan. 

Received: 11 June 2013 Accepted: 19 June 2013 
Published: 27 June 2013 



Ikuta et al. BMC Evolutionary Biology 201 3, 13:1 29 
http://www.biomedcentral.eom/1 471 -21 48/1 3/1 29 



Page 9 of 10 



References 

1. Holland PW: Beyond the Hox: how widespread is homeobox gene 
clustering? J Anat 2001, 199:13-23. 

2. Ferrier DE, Holland PW: Ancient origin of the Hox gene cluster. Nat Rev 
Genet 2001,2:33-38. 

3. Garcia-Fernandez J: The genesis and evolution of homeobox gene 
clusters. Nat Rev Genet 2005, 6:881-892. 

4. Hui JH, Holland PW, Ferrier DE: Do cnidarians have a ParaHox cluster? 
Analysis of synteny around a Nematostella homeobox gene cluster. 
Evol Dev 2008, 10:725-730. 

5. Mendivil Ramos 0, Barker D, Ferrier DE: Ghost loci imply Hox and ParaHox 
existence in the last common ancestor of animals. Curr Biol 2012, 
22:1951-1956. 

6. Ryan JF, Mazza ME, Pang K, Matus DQ, Baxevanis AD, Martindale MQ, 
Finnerty JR: Pre-bilaterian origins of the Hox cluster and the Hox code: 
evidence from the sea anemone, Nematostella vectensis. PloS one 2007, 
2(1):e153. 

7. Ferrier DE, Minguillon C: Evolution of the Hox/ParaHox gene clusters. 

IntJDev Biol 2003, 47:605-611. 

8. McGinnis W, Krumlauf R: Homeobox genes and axial patterning. Cell 1992, 
68:283-302. 

9. Izpisua-Belmonte JC, Falkenstein H, Dolle P, Renucci A, Duboule D: Murine 
genes related to the Drosophila AbdB homeotic genes are sequentially 
expressed during development of the posterior part of the body. EMBO J 
1991, 10:2279-2289. 

10. Duboule D, Morata G: Colinearity and functional hierarchy among genes 
of the homeotic complexes. Trends Genet 1994, 10:358-364. 

1 1 . Brooke NM, Garcia-Fernandez J, Holland PW: The ParaHox gene cluster is 
an evolutionary sister of the Hox gene cluster. Nature 1998, 392:920-922. 

12. Ferrier DE, Dewar K, Cook A, Chang JL, Hill-Force A, Amemiya C: The 
chordate ParaHox cluster. Curr Biol 2005, 15:R820-R822. 

13. Furlong RF, Mulley JF: ParaHox cluster evolution-hagfish and beyond. 
Zoolog Sci 2008, 25:955-959. 

14. Weiss JB, Von Ohlen T, Mellerick DM, Dressier G, Doe CQ, Scott MP: 
Dorsoventral patterning in the Drosophila central nervous system: the 
intermediate neuroblasts defective homeobox gene specifies 
intermediate column identity. Genes Dev 1998, 12:3591-3602. 

15. Macdonald PM, Struhl G: A molecular gradient in early Drosophila 
embryos and its role in specifying the body pattern. Nature 1986, 
324:537-545. 

16. Ruvkun G, Hobert O: The taxonomy of developmental control in 
Caenorhabditis elegans. Science 1 998, 282:2033-2041 . 

17. Hui JH, Raible F, Korchagina N, Dray N, Samain S, Magdelenat G, Jubin C, 
Segurens B, Balavoine G, Arendt D, et al: Features of the ancestral 
bilaterian inferred from Platynereis dumerilii ParaHox genes. BMC Biol 
2009, 7:43. 

18. Mulley JF, Chiu CH, Holland PW: Breakup of a homeobox cluster 
after genome duplication in teleosts. Proc Natl Acad Sci USA 2006, 
103:10369-10372. 

19. Ferrier DE, Holland PW: Ciona intestinalis ParaHox genes: evolution of 
Hox/ParaHox cluster integrity, developmental mode, and temporal 
colinearity. Mol Phylogenet Evol 2002, 24:41 2-41 7. 

20. Arnone Ml, Rizzo F, Annunciata R, Cameron RA, Peterson KJ, Martinez P: 
Genetic organization and embryonic expression of the ParaHox genes in 
the sea urchin S. purpuratus: insights into the relationship between 
clustering and colinearity. Dev Biol 2006, 300:63-73. 

21. Annunziata R, Martinez P, Arnone Ml: Intact cluster and conserved 
expression of ParaHox genes in the asteroid Patiria miniata. BMC Biol 
2013, 11:68. 

22. Osborne PW, Benoit G, Laudet V, Schubert M, Ferrier DE: Differential 
regulation of ParaHox genes by retinoic acid in the invertebrate 
chordate amphioxus {Branchiostoma floridae). Dev Biol 2009, 
327:252-262. 

23. Hudson C, Lemaire P: Induction of anterior neural fates in the ascidian 
Ciona intestinalis. Mech Dev 2001, 100:189-203. 

24. Corrado M, Aniello F, Fucci L, Bran no M: Ci-IPFl, the pancreatic 
homeodomain transcription factor, is expressed in neural cells of Ciona 
intestinalis larva. Mech Dev 2001, 102:271-274. 

25. Kusakabe T, Yoshida R, Kawakami I, Kusakabe R, Mochizuki Y, Yamada L, 
Shin-i T, Kohara Y, Satoh N, Tsuda M, et al: Gene expression profiles in 
tadpole larvae of Ciona intestinalis. Dev Biol 2002, 242:188-203. 



26. Imai KS, Hino K, Yagi K, Satoh N, Satou Y: Gene expression profiles of 
transcription factors and signaling molecules in the ascidian embryo: 
towards a comprehensive understanding of gene networks. Development 
2004, 131:4047-4058. 

27. Sherwood Rl, Chen TY, Melton DA: Transcriptional dynamics of 
endodermal organ formation. Dev Dyn 2009, 238:29-42. 

28. Li H, Arber S, Jessell TM, Edlund H: Selective agenesis of the dorsal 
pancreas in mice lacking homeobox gene Hlxb9. Nat Genet 1999, 
23:67-70. 

29. Valerius MT, Li H, Stock JL, Weinstein M, Kaur S, Singh G, Potter SS: Gsh-1: a 
novel murine homeobox gene expressed in the central nervous system. 
Dev Dyn 1995, 203:337-351. 

30. Frobius AC, Matus DQ, Seaver EC: Genomic organization and expression 
demonstrate spatial and temporal Hox gene colinearity in the 
lophotrochozoan Capitella sp. I. PloS One 2008, 3:e4004. 

31. Frobius AC, Seaver EC: ParaHox gene expression in the polychaete 
annelid Capitella sp. I. Dev Genes Evol 2006, 216:81-88. 

32. Samadi L, Steiner G: Conservation of ParaHox genes' function in 
patterning of the digestive tract of the marine gastropod Gibbula varia. 
BMC Dev Biol 2010, 10:74. 

33. Ikuta T: Evolution of invertebrate deuterostomes and Hox/ParaHox 
genes. Genomics Proteomics Bioinformatics 201 1, 9:77-96. 

34. Rottinger E, Lowe CJ: Evolutionary crossroads in developmental biology: 
hemichordates. Development 2012, 139:2463-2475. 

35. Brown FD, Prendergast A, Swalla BJ: Man Is but a Worm: Chordate Origins. 
Genesis 2008, 46:605-613. 

36. Caron JB, Morris SC, Cameron CB: Tubicolous enteropneusts from the 
Cambrian period. Nature 2013, 495:503-506. 

37. Cameron CB, Garey JR, Swalla BJ: Evolution of the chordate body plan: 
new insights from phylogenetic analyses of deuterostome phyla. Proc 
Natl Acad Sci USA 2000, 97:4469-4474. 

38. Cannon JT, Rychel AL, Eccleston H, Halanych KM, Swalla BJ: Molecular 
phylogeny of hemichordata, with updated status of deep-sea 
enteropneusts. Mol Phylogenet Evol 2009, 52:1 7-24. 

39. Peterson KJ, Su YH, Arnone Ml, Swalla BJ, King B: microRNAs support the 
monophyly of enteropneust hemichordates. J Exp Zool B 2013. 
doi:10.1002/jez.b.22510. 

40. Byrne M, Nakajima Y, Chee FC, Burke RD: Apical organs in echinoderm 
larvae: insights into larval evolution in the Ambulacraria. Evol Dev 2007, 
9:432-445. 

41. Harada Y, Shoguchi E, Taguchi S, Okai N, Humphreys T, Tagawa K, Satoh N: 
Conserved expression pattern of BMP-2/4 in hemichordate acorn worm 
and echinoderm sea cucumber embryos. Zoolog Sci 2002, 1 9:1 1 1 3-1 1 21 . 

42. Rottinger E, Martindale MQ: Ventralization of an indirect developing 
hemichordate by NiCI suggests a conserved mechanism of dorso-ventral 
(D/V) patterning in Ambulacraria (hemichordates and echinoderms). 
Dev Biol 2011, 354:173-190. 

43. Freeman R, Ikuta T, Wu M, Koyanagi R, Kawashima T, Tagawa K, 
Humphreys T, Fang GC, Fujiyama A, Saiga H, et al: Identical genomic 
organization of two hemichordate hox clusters. Curr Biol 2012, 
22:2053-2058. 

44. Peterson KJ: Isolation of Hox and Parahox genes in the hemichordate 
Ptychodera flava and the evolution of deuterostome Hox genes. 

Mol Phylogenet Evol 2004, 31:1 208-1 21 5. 

45. Ikuta T, Miyamoto N, Saito Y, Wada H, Satoh N, Saiga H: Ambulacrarian 
prototypical Hox and ParaHox gene complements of the indirect- 
developing hemichordate Balanoglossus simodensis. Dev Genes Evol 2009, 

219:383-389. 

46. Simakov O, Marletaz F, Cho SJ, Edsinger-Gonzales E, Havlak P, Hellsten U, 
Kuo DH, Larsson T, Lv J, Arendt D, et al: Insights into bilaterian evolution 
from three spiralian genomes. Nature 2013, 493:526-531. 

47. Zhang G, Fang X, Guo X, Li L, Luo R, Xu F, Yang P, Zhang L, Wang X, Qi H, 
et al: The oyster genome reveals stress adaptation and complexity of 
shell formation. Nature 2012, 490:49-54. 

48. Takeuchi T, Kawashima T, Koyanagi R, Gyoja F, Tanaka M, Ikuta T, 
Shoguchi E, Fujiwara M, Shinzato C, Hisata K, et al: Draft genome of the 
pearl oyster Pinctada fucata: a platform for understanding bivalve 
biology. DNA Res 2012, 19:117-130. 

49. Fujiyama A, Watanabe H, Toyoda A, Taylor TD, Itoh T, Tsai SF, Park HS, 
Yaspo ML, Lehrach H, Chen Z, et al: Construction and analysis of a 
human-chimpanzee comparative clone map. Science 2002, 295:131-134. 



Ikuta et al. BMC Evolutionary Biology 201 3, 13:1 29 
http://www.biomedcentral.eom/1 471 -21 48/1 3/1 29 



Page 10 of 10 



50. Tagawa K, Nishino A, Humphreys T, Satoh N: The spawning and early 
development of the Hawaiian acorn worm (Hemichordate), Ptychodera 

flava.Zoolog Sci 1998, 15:85-91. 

51. Lowe CJ, Tagawa K, Humphreys T, Kirschner M, Gerhart J: Hemichordate 
embryos: procurement, culture, and basic methods. Methods Cell Biol 
2004, 74:171-194. 

52. Winchell CJ, Sullivan J, Cameron CB, Swalla BJ, Mallatt J: Evaluating 
hypotheses of deuterostome phylogeny and chordate evolution with 
new LSU and SSU ribosomal DNA data. Mol Biol Evol 2002, 19:762-776. 

53. Chen JH, Luo YJ, Su YH: The dynamic gene expression patterns of 
transcription factors constituting the sea urchin aboral ectoderm gene 
regulatory network. Dev Dyn 201 1, 240:250-260. 



doi:1 0.1 1 86/1 471 -21 48-1 3-1 29 

Cite this article as: Ikuta et al.: Identification of an intact ParaHox cluster 
with temporal colinearity but altered spatial colinearity in the 
hemichordate Ptychodera flava. BMC Evolutionary Biology 2013 13:129. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at (^\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ ™omea centra I 



